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Abstract. We investigate the RPA normal-mode coordinates in the pairing- 
plus-quadrupole model, with an eye on simplifying the application of large 
amplitude collective motion techniques. At the Hartree-Bogoliubov minimum, 
Os . the RPA modes are exactly the cranking operators of the collective coordinate 

approach. We examine the possibility of representing the self-consistent cranking 
operator by linear combinations of a limited number of one-body operators. We 
^ ' study the Sm nuclei as an example, and find that such representations exist in 

I/"") , terms of operators that are state-dependent in a characteristic manner. 
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Q 1 PACS numbers: 21.60.-n, 21.60.Ev, 21.60.Jz 

Q\ , The selection of proper collective variables is an important problem in the study 

OS ' of large amplitude collective motion. In the usual Constrained Hartree-Fock (CHF) or 

Hartree-Fock-Bogoliubov (CHFB) calculations, the collective subspaces are generated 
by a small number of one-body constraint (also called cranking) operators which are 
most commonly taken to be of the multipole form (t l Ylk)- In realistic calculations 
of processes such as fission, the number of coordinates to describe the full nuclear 
dynamics can easily become larger than can be dealt with in satisfactory manner, and 
a method to determine the optimal combination needs to be devised. Even assuming 
that such a method exists, there is no a priori reason to limit oneself to multipole 
operators, and the cranking operators should be determined by the nuclear collective 
dynamics itself, from the set of all one-body operators. This is clearly a difficult task, 
and one would like to be able to select a small group of operators, and find the optimal 
combination of these operators at each point of the collective surface. 

In our past work, we have investigated a theory of adiabatic large amplitude 
collective motion as a method to generate self-consistent collective subspaces (see 
reference and references therein). The key ingredient of the method is the self- 
consistent determination of the constraint operator, and as such it may provide an 
answer to the selection question discussed above. Using the local harmonic version 
(LHA) of the theory §, we have recently embarked on a study of the properties of large 
amplitude collective motion in systems with pairing. We have dealt with two simple 
models: a semi-microscopic model of nucleons interacting through a pairing force, 
coupled to a single harmonic variable M and a fully microscopic 0(4) model which may 
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be regarded as a simplified version of the pairing-plus-quadrupole (P+Q) Hamiltonian 
It has turned out that the self-consistent collective coordinate obtained by the 
LHA accounts quite well for the exact dynamics of these models. We have also shown 
that the CHFB calculations using the mass-quadrupole operator as the constraint 
operator can result in incorrect results || . It is not immediately obvious that we can 
extrapolate the conclusions reached in these models to fully realistic nuclear problems. 
Within the time-dependent Hartree-Bogoliubov (TDHB) approximation, the dynamics 
of these models is described by a modest number of degrees of freedom (4 ~ 12). For 
realistic problems in heavy nuclei, on the other hand, we need to deal with millions of 
degrees of freedom! In this letter, we report the first attempts to study such realistic 
nuclear problems. 

In the LHA, the collective path (or collective surface for more than one coordinate) 
is determined by solving the CHFB problem with a cranking operator which is self- 
consistently determined by the local RPA. Since this procedure requires us to solve the 
RPA at each point on the collective path, it will be very useful if the RPA eigenvectors 
can be approximated by taking linear combinations of a small number of one-body 
operators. This means that we can restrict the RPA diagonalisation to this small 
space, rather than deal with the full millions-by-millions RPA matrix. Such a scenario 
has been examined in reference Q] and in reference ^ for the HF problem in 28 Si, and 
a strong radial and spin dependence of self-consistent cranking operators was found. 
However, since the model space had only six degrees of freedom and the pairing degrees 
of freedom were neglected, the results can not be directly generalised to heavy nuclei. 
Here we perform a similar analysis for the P+Q model. Since the model is known 
to be able to realistically describe collective phenomena involving both pairing and 
quadrupole degrees of freedom || fjj, we expect that the same choice for the set of 
one-body operators should work for other realistic Hamiltonians as well. 

Baranger and Kumar analysed in great detail the collective motion in the P+Q 
model assuming that the collective variables are the mass quadrupole operators JtJ. 
Thus, they reduced the large number of two-quasiparticle (2qp) degrees of freedom (of 
the order of a thousand) into only two "collective" coordinates, (5 and 7. However, our 
previous study of the 0(4) model || suggests that even for such simple Hamiltonians 
the self-consistent collective coordinate is not as trivial as it seems to be. In this 
letter we report the first result for the P+Q model and show that the normal-mode 
coordinate of the random-phase approximation (RPA) is quite different from the mass 
quadrupole operator. This is especially true when the system is deformed. 

For the P+Q model, the Hartree-Bogoliubov (HB) ground state can be specified 
by 6 parameters (A„, A p , A„, X p , (3 and 7). The TDHB equations can be shown to be 
equivalent to Hamilton's equations of motion [|| . The underlying classical Hamiltonian 
can be written, up to second order in coordinates and momenta, as 

H = (9\H - nN\9) ^E + -B^irp + \c a& Cf , (1) 

in terms of the canonical variables (£, n) B. Here each of the indices (a,/3, •••) 
indicates a pair of 2qp indices (ij,kl, ■ ■ •). We adopt the standard tensor notation 
where a repeated upper and lower index indicates a summation. The mass and 
curvature parameters are 

B af) = E a S aP - 2 ]T XpS^S^ , C aP = E a 6 a0 - 2 £ x P R { a p) R? ■ (2) 
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Here R 1 -^ and represent the hermitian and anti-hermitean components, 

respectively, of the pairing and quadrupole operators, 

R^ = ^2R^((^^) a +h.c.) , SW =J2s { a p) ((aV) Q -h.c.) , (3) 

a a 

where the scattering terms (a^o) are omitted. Following reference Q, we multiply the 
quadrupole operators by a factor a 2 , with a n = {2N/A) 1 / Z and a p = (2Z ' j 'A) 1 / 3, , and 
also reduce the quadrupole matrix elements between the states of the upper shell by 
a factor Q = (A/l + §)/(A/" + |), where M is the oscillator quantum number operator 
and jVl is the number of quanta in the lower shell. Thus, the modified quadrupole 
operators are defined as Q 2K = {Q2k)u + (Q2k) p , with (Q 2 k)u = c? n (,(r 2 Y 2K ) n 
and (Q 2 k) p = ce 2 ((r 2 Y 2 K) p (which we shall refer to as "the quadrupole operators"). 
The RPA equation is solved by diagonalisation of the Hamiltonian (|l|). This can be 
implemented by a (linear) transformation to normal coordinates (q^), 

q» = f%e, r = (4) 

where f^ a and g ^ are just coefficients of linear transformation. These can be also 
identified as dq^ 1 /d^ a and d^ a jdq 11 for the general point transformation, = /^(t;) 
and £ a = g a {q), when we study large excursions from equilibrium Q. In the case of a 
single collective coordinate, the path is determined by solving the CHB equation with 
the cranking operator q 1 — / 1 (^)- A similar set of equations has been also obtained 
by other adiabatic time-dependent mean-field theories . 
The solution of the RPA equation 

c a7 ir*»/£ = m 2 f% , (5) 

involves the diagonalisation of the RPA matrix C Q7 -B 7 ^ whose dimension is equal 
to the number of active 2qp degrees of freedom. For separable forces, this can 
be simplified by solving a dispersion relation. In general, however, the RPA 
diagonalisation requires extensive computational resources. Now let us approximate 
an eigenvector using a certain set of one-body operators {OW}: 

£« = £c,o« , (6) 

i 

where O^a indicate the 2qp matrix elements of operator OW as in equation (S). Then, 
instead of the full RPA equation (||), we obtain a projected RPA equation 

M ij C? = (n n ) 2 N ij CJ , (7) 

which determines the coefficients C™, with 

M ij = O^B^C^B^O^ , N ij = O^B^of . (8) 

The dimension of the matrices M y and A™ is equal to the number of one-body 
operators {OW}. Therefore, if we can approximate the RPA eigenvectors by using a 
small number of operators, it will significantly reduce the computational task. 

A criterion for good projection may be given by the closeness of the projected 
RPA frequencies Ct to the real RPA frequencies il. Another criterion is the smallncss 
of the quantity 5, 

5 = !, a B a(i (f tP - f J3 ) . (9) 

If we normalise / p as f^B^fg = , we find < 5 < 1 with 6 = corresponding 
to the exact projection and 5 = 1 to the case where / is orthogonal to /. 
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Figure 1. Calculated excitation energies of quadrupole, j3 and 7 vibrations for 
even-even Sm isotopes. Note that the ground states of 146,148 Sm are spherical. 
The closed circles indicate the RPA results while the open symbols are the results 
of projected RPA calculations. See the main text for the difference between (1), 
(2) and (3). 

We have performed calculations for several heavy isotopes. Here we report the 
numerical results for even-even Sm isotopes (A=146~154). The form of the P+Q 
model is that discussed in the second and third of the series of papers by Baranger 
and Kumar H|. The model space and the parameters, such as the spherical single- 
particle energies, the pairing and quadrupole force strengths, are taken from table 1 
in the third paper. The equilibrium parameters (/3, 7, A, A) are found to agree with 
table 2 of the same paper. The ground states of 146 ' 148 Sm are spherical (f3 = 0) and 
the others have prolate shapes (/3 > 0,7 = 0). 

Figure [l] shows the excitation energies (RPA frequencies) of (3 and 7 vibrations, 
obtained by the RPA and projected RPA calculations. For the projected RPA 
calculations, we have adopted three different sets of one-body operators. The first 
and simplest choice is to use the operators appearing in the separable forces, the 
pairing and quadrupole operators, P T , PJ, (Q2k)t (t = n,p). This choice is denoted 
as (1) in the figure. In this case the projected RPA matrices of equation (||) are two 
dimensional for spherical nuclei and for 7 vibrations, while they are six dimensional 
for the (3 vibrations. The calculated frequencies are 7 ~ 8 MeV which are 5^6 
MeV larger than the corresponding RPA frequencies. In the second set, labelled as 
(2), we increase the number of operators. We keep both the pairing and quadrupole 
operators, but include two additional quadrupole operators with "monopole" radial 
dependence, (r Y 2 K)T- We also include the hexadecapole operators, (V 4 >4A')t, and 

the rank-2 spin-dependent operators, ([r !^ x s ]k)t' (I 7 * 2 -^ x s ]a')t- As ^ ar as * ne 
frequencies are concerned, we can see some improvement over the case (1) for spherical 
and the 7 vibrations in deformed nuclei, though they are still much higher than the 
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Table 1. Calculated values of 5, equation (g), for the projected RPA solutions 
for Sm isotopes. The columns (1), (2), (3), (1-a) and (1-b) represent the different 
projections (see text). For the spherical nuclei ( 146,148 Sm), there is no distinction 
between and 7 vibrations. 



A (3 vibration 7 vibration 





(1) 


(1-a) 


(1-b) 


(2) 


(3) 


(1) 


(1-a) 


(1-b) 


(2) 


(3) 


146 


0.271 


0.132 


0.225 


0.421 


0.009 












148 


0.243 


0.131 


0.184 


0.314 


0.0003 












150 


0.602 


0.499 


0.519 


0.632 


0.026 


0.610 


0.342 


0.507 


0.685 


0.092 


152 


0.497 


0.346 


0.433 


0.526 


0.020 


0.616 


0.279 


0.472 


0.691 


0.081 


154 


0.513 


0.117 


0.437 


0.534 


0.002 


0.636 


0.208 


0.426 


0.679 


0.052 



real RPA frequencies. For the (3 vibrations, the inclusion of the additional rank-2 (and 
higher rank) operators seems not so important. Actually we see that the (3 vibrations 
are found to have a significant amount of monopole components. For the last set, 
denoted as (3), we adopt the same operators as (1) but each 2qp matrix element is 
weighted with a factor (-E^qp) -2 - This means that we employ a set of state- dependent 
one-body operators {OW} defined by 

5 ^E?|y( atat )«+ h - c -- ( 10 ) 

The result of this projection is now almost identical to that of the full RPA. 

In table |l], the quality of projection 5, equation ([}]), is listed. In the cases (1) 
and (2), where the RPA vectors are projected on the elementary operators, 5 > 0.25 
for 146,148 Sm and S > 0.5 for the others. Therefore, roughly speaking, the one-body 
operators possess at most 75% (50%) of overlap with the real eigenvectors in spherical 
(deformed) nuclei. On the other hand, the projection (3) exhausts more than 90% of 
real eigenvectors even for the worst case. At first sight it may look strange that 6 is 
larger for (2) than for (1), while the energy for (2) is lower. This is due to the fact 
that case (2) is dominated by certain neutron components. Since the relevant neutron 
2qp energies are lower than those of protons, this proton- neutron asymmetry leads to 
a decrease in the frequency f2 and at the same time an increase in 5. This is also a 
reflection of the poor quality of the approximation. 

The figure and table indicate that it is very difficult to obtain sensible results by 
using elementary one-body operators (i.e., of the form (1) or (2)). This is mainly due 
to the fact that the RPA eigenvectors, when being projected onto elementary one- 
body operators, have unrealistically large amplitudes for high-lying 2qp components. 
In order to demonstrate this, we introduce a cut-off energy A cut for the 2qp matrix 
elements, i.e., O^a = for E a > A cut . We then perform the projected RPA calculation 
with the truncated one-body operators from set (1). The resulting values 5 are listed 
in table for A cut = 5 MeV (1-a) and for 10 MeV (1-b). We see that the major 
contributions to the RPA modes come from the 2qp components with -E^qp < 5 MeV. 
We thus conclude that the superiority of the projection (3) simply comes from its 
being capable of suppressing the unnecessary high-energy components by the factor 
(E^qp) -2 - This suppression factor is not arbitrary, but can be derived from the 
following simple argument. If we have a single-mode separable force H = —{\/2)xRR 



Letter to the Editor 



(> 



(assuming no coupling among different modes), we can analytically determine the RPA 
eigenvectors, oc R a /((E a ) 2 — (f^) 2 ). In the limit that <C -E^qp, the projection 
on f a oc R a /(E a ) 2 gives the exact answer. 

Let us examine the projection (3) in more detail. For spherical Sm nuclei, the RPA 
eigenvector is of isoscalar character and can be approximated as / ss (Q2) n + {Q2)p 
where the tilde indicates that the matrix elements include the suppression factor as in 
equation ([lO]). For deformed nuclei, where the collectivity of the vibrational states is 
smaller than for spherical nuclei and the pairing modes can mix with the quadrupolc 
ones, the situation is more complex. Taking 154 Sm as an example, the eigenvectors of 
P and 7 vibrations are 

j/3-vib = ( g 2())n + o.9i (Q 20 ) p - .48F„ - O.UP p + 0.085F,! - O.FLPj , (11) 
fy-va, = { Q 22)n + .87(Q 22 ) P . (12) 

For the (3 vibration, we find a significant mixing with the monopole pairing modes. 

In conclusion, we have examined the possibility of expressing the self-consistent 
cranking operator in terms of a limited set of one-body operators. It seems very 
difficult to approximate the normal-mode vectors with use of elementary one-body 
operators. This difficulty disappears, however, when we use a small number of 
state- dependent one-body operators. This may reflect the importance of the self- 
consistent determination of the collective coordinates for large amplitude collective 
motion, because the coordinates now have a strong state-dependence as well. The 
structure of the self-consistent cranking operators is clearly changing when we move 
from spherical to deformed nuclei. For the study of large amplitude collective motion 
in heavy nuclei for which the diagonalisation of the RPA matrix becomes too time- 
consuming, the results of this paper may give a hint for a correct choice of a state- 
dependent basis of operators. The choice of a limited set of (state dependent) basis 
operators provides a practical way to solve the LHA through the projection. With the 
self-consistent cranking operators, the LHA should provide a significant improvement 
over the conventional CHFB calculation based on fixed cranking operators. Clearly 
we have not discussed the structure of the self-consistent cranking operator away from 
the minimum point. This will be the subject of a future publication. 
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